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The  study  of  resonant  tunneling  structures  (RTS)  has  become  increasingly  important  be¬ 
cause  of  possible  device  applications  and  also  because  of  the  basic  physics  they  involve. 
Within  the  static  picture  we  have  calculated  the  dwell  time,  i.e.,  the  average  time  an 
electron  spends  inside  an  RTS.  The  actual  time  involved  in  the  tunneling  process  can 
be  described  by  the  solution  of  the  time-dependent  effective  mass  Schodinger  equation 
(TDEMSE),  and  a  code  based  on  the  fast  Fourier  transform  method  is  developed  to  solve 
the  TDEMSE  for  a  potential  profile  such  as  RTS.  The  physical  properties  of  an  RTS  in¬ 
vestigated  are  the  build-up  time,  which  is  the  time  taken  by  an  electron  to  accumulate 
probability  inside  the  well,  and  the  escape  time. 


1.  Introduction 

How  long  it  takes  for  a  moving  particle  to  penetrate  a  potential  barrier  is  a  question 
being  asked  often  since  the  discovery  of  alpha  particles.  Recently,  interest  in  investigating 
the  time  scale  involved  in  the  tunneling  process  was  revived  after  the  work  of  Tsu  and 
Esaki.1  where  they  explained  the  negative  differential  resistance  in  the  I —  V  characteristic 
of  a  compound  semiconductor  diode  through  the  resonant  tunneling  process.  Since  then, 
several  statical  and  dynamical  physical  properties  of  resonant  tunneling  structures  (a  single 
well  formed  by  two  barriers  of  finite  widths  and  heights.2  henceforth  referred  to  as  RTS) 
have  been  the  subject  of  investigation.  A  vast  literature3'4  is  available  by  now  on  the  field. 
Our  main  interest  here  is  to  answer  the  question  above,  although  the  time  scale  involved  in 
the  resonant  tunneling  process  is  the  most  controversial  topic  still  today.  We  have  studied 
through  a  static  picture  the  most  accepted  dynamical  quantity,  the  dwell-time,  which  is 
the  average  time  spent  by  an  electron  inside  an  RTS.  A  dynamical  aspect  of  the  problem 
is  also  addressed  by  the  solution  of  the  one- dimensional  time  dependent  effective  mass 
Schrodinger  equation  (TDEMSE),  and  the  properties  investigated  from  the  solution  are 
the  build-up  time,  which  is  the  time  taken  by  an  electron  to  accumulate  probability  inside 
the  well,  and  the  escape  time,  which  is  on  a  time  scale  associated  with  the  resonance  width 
through  the  uncertainty  principle. 

2.  Static  Picture  and  Dwell-Time 

As  mentioned  above,  the  dwell  time  is  an  average  time  during  which  the  probability 
to  find  the  electron  inside  an  RTS  is  finite.  Hence,  the  dwell  time  is  the  total  probability 
of  the  electron  per  unit  incident  flux.0  The  stationary-state  properties  of  an  RTS  in  the 
effective  mass  approximation  are  obtained  by  solving  the  time-independent  Schrodinger 
equation  for  the  envelop  function  'l'(r)  along  the  growth  direction  i  :6'7 

+  lV(x)  ~  =  °  i 


(1) 


.  2  . 


where  V(x )  is  the  RTS  potential  profile 

(  V i ,  if  0  <  x  <  ai  , 

V(x)  =  <  0,  if  a\  <  x  <  aj  +  d  , 

t  V2,  if  dj  4"  d  <  x  <  ai  +  d.  +  02  , 

V{  and  di  (i— 1,2)  are  the  heights  and  widths  of  the  barriers,  d  is  the  width  of  the  well, 
m  is  effective  mass,  and  2a  +  6  =  —1  with  a  and  6  as  constants.  These  constants  must 
appear  in  the  form  given  in  order  to  keep  the  Hamiltonian  Hermitian.  The  kinetic  energy 
operator  of  Eq.  (1)  dictates  that  ma'I'(x)  and  ma+b^j—  must  be  continuous  across  the 
interface,  implying  the  physical  result  that  the  current  density  j  «  ^f-77  =  ma'I,*m0+k^- 
be  continuous.  However,  in  general,  the  charge  density  p  oc  need  not  be  continuous 
across  an  interface.  For  the  special  case  of  a  =  0  and  b  =  —  1  one  obtains,  in  addition,  the 
continuity  of  charge  density.  So,  this  is  a  single-parameter  b  (a  =  —  problem.  The 
resonance  states  clearly  depend  on  this  parameter  6,  and  one  can  not.  a  priori ,  prefer  one 
value  of  b  over  another. 

We  have  studied  the  dependence  of  the  dwell  time  and  width  of  the  resonance  on  b. 
The  dwell  time  tq  over  the  region  of  0  to  xj  of  the  structure  is  defined4  as  the  integrated 
probability  density  of  the  electron  per  unit  incidence  flux, 

T71  ^ 

TD=  hk  Jo  dx  ^  I2  • 

where  k  =  yj Figures  1  and  2.  and  3  and  4  show  the  dwell  time  as  a  function  of  x  for 
three  different  6  parameter  values,  -2,  0  &  2,  and  width  of  the  resonance  as  a  function  of 
b  for  the  first  and  second  resonance  states,  respectively.  It  is  clear  from  Figs.  1  and  3  that 
the  dwell  rime  varies  considerably  as  one  changes  6,  but  a  reverse  pattern  is  found  (Figs. 
2  and  4)  in  the  change  of  resonance  widths  of  the  first  and  second  resonances. 

We  have  thus  shown  that  the  parameter  associated  with  the  boundary  conditions 
has  a  profound  effect  on  the  characteristics  of  the  system  such  as  resonance  energy  and 
width,  and  dwell  time.  It  has  been  shown8  in  a  different  context  that  for  two  semi¬ 
infinite  heterostructures,  the  boundary  conditions  at  the  interface  not  only  involve  the 
effective  masses  but  also  certain  other  parameters  which  are  microscopic  in  origin,  having 
no  macroscopic  analogs.  We  therefore  believe  that  the  arbitrariness  in  the  choice  of  the 
parameter  b  can  be  fixed  through  microscopic  calculations. 

2.  Dynamic  Picture  and  Build-up  and  Escape  Times 

For  b  =0  and  constant  mass  throughout  the  RTS,  Eq.  (1)  with  E  =  ih^  is  solve 
numerically  with  the  fast  Fourier  transform  method.9  When  using  this  method,  it  has  been 
found  convenient  to  model  a  discontinuous  potential,  such  as  an  RTS.  as  a  combination  of 
Fermi  distribution  functions  of  the  form 

V(x)  _  _ 1 _ 1 _ 

Vq  exp[(x  -  x0  -  cj2)/S]  +  1  exp[(x  -  x0  +  c/2)/£]  +  1 

_ 1 _ t  .  ^ 

exp[(x  -  x0  -  d/2)/«5]  +  1  exp[(x  -  x0  +  d/2)/ 6}  4- 1 
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Fig.  1.  Dwell-time  as  a  function  of  posi¬ 
tion  x  for  the  first  resonance  state.  The 
solid,  dashed  and  dot-dashed  axe  for 
b  —  *2,  0  and  2,  respectively.  The  other 
parameters  are:  ax  =  <z2=50  A,  d=  100 


Fig.  2.  Width  of  the  first  resonance  state 
as  a  function  of  b. 
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Fig.  3.  Same  as  in  Fig.  1  but  for  the 
second  resonance  state. 
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Fig.  4.  Width  of  the  second  resonance 
state  as  a  function  of  6. 
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where,  Vo  =  V\  =  Vi  is  the  barrier  height,  c  =  a\+d  +  ai  is  the  total  length  of  the  RTS,  xo 
is  the  position  of  the  middle  of  the  RTS,  and  8  is  a  smoothing  parameter.  We  have  carried 
out  a  single-step  calculation  using  1024  grid  points.  20  out  of  which  are  used  to  discretize 
the  RTS,  with  10  out  of  this  20  to  discretize  the  well  in  the  RTS.  Here  8  is  taken  to  be  2 
A,  and  a  box  of  length  10240  A  corresponding  to  a  grid  spacing  of  10  A  has  been  chosen 
for  this  calculation.  These  parameters  are  found  appropriate  to  demonstrate  the  dynamics 
of  an  electron  in  a  box  which  has  a  symmetric  RTS  located  at  its  center. 

To  study  the  transport  properties  of  an  electron,  we  have  chosen  its  initial  wavefunc- 
tion  to  be  represented  by  a  Gaussian  wavepacket, 

=  0)  =  N exp[— — 1  — -  +  ik0x]  ,  (4) 

4  o* 

where  xp  is  the  center  of  the  wavepacket  (one  fourth  of  the  box  length  from  the  left  wall), 
a  is  the  width  of  the  wavepacket.  ko  is  the  mean  wave  vector  (A k  =  £),  and  N  is  a 
normalization  factor.  The  wavefunction  '$(x>t)  has  been  calculated  for  different  values  of 
t. 

A  symmetric  RTS  with  c  =  200  A,  d  —  100  A  and  Vo  =  200  meV  is  analyzed.  There 
are  two  resonance  states  at  30.87  ±0.174  meV  and  117.4  ±2.98  meV  of  the  RTS  considered 
here.  Figures  5  and  6  display  the  integrated  probability  density  P(t)  =  dx|^(x,  t)\2 

inside  the  well,  which  is  part  of  the  wavefunction  trapped  in  the  well  as  a  function  of 
time  for  the  first  and  second  resonance  states,  respectively.  The  transmission  probability 
T(t)  =  f^+c/2  dx|'I,(x,i)|2  is  also  shown  in  Figs.  5  and  6  by  dotted  lines.  It  is  clear 
that  P(t)  increases  up  to  a  certain  value  Pmax,  and  the  time  taken  to  reach  Pmax  from  a 
reference  point  is  called  the  buildup-time  rj.  The  probability  P(t)  decays  exponentially 
after  the  build-up  time,  and  an  exponential  fit  to  the  decay  part  of  P(t)  versus  time, 
P(t)  =  Pmax  exp(— Qt),  yields  0=0.264  and  4.5  ps~l  and  hence,  the  width  of  the  state 
Aj Er  =  0.174  and  2.98  meV  for  the  first  and  second  resonance  states,  respectively.  Since 
the  first  resonance  state  is  sharp  in  energy,  it  takes  more  time  to  decay,  during  which  the 
reflected  and  transmitted  parts  of  the  wavepacket  return  back  from  the  box  wall.  To  avoid 
interferences  with  the  returned  reflected  and  transmitted  parts  of  the  wavepacket,  instead 
of  allowing  P(t)  to  decay  further,  a  fitted  curve  is  shown  by  the  dot-dashed  line  in  Fig.  5. 

In  Figure  7,  we  show  (a)  the  maximum  probability  density  Pma x  trapped  in  the 
well  in  the  upper  panel  and  (b)  the  transmission  probability  along  with  the  plane  wave 
transmission  coefficient  (solid  line)  in  the  lower  panel  as  a  function  of  the  electron  energy 
in  the  neighborhood  of  the  second  resonance  state  for  a  —  400  (solid  circles),  500  (open 
circles)  and  600  (solid  triangles)  A.  The  lower  panel  also  displays  the  exact  transmission 
coefficient  T9  for  a  wavepacket10  with  cr=  400  (dashed  line),  500  (dotted  line)  and  600 
(dot-dashed  line)  A.  The  resonance  energy  of  T{t)  in  the  lower  panel  of  Fig.  7  is  slightly 
shifted  from  the  T9  because  of  the  reason  that  the  Fermi  distributed  structure  as  in  Eq. 
(3)  is  not  a  exact  representation  of  the  RTS  used  to  calculate  the  plane  wave  transmission 
coefficient  and  T<T.  The  transmission  probability  T(t)  associated  with  the  dynamics  and 
exact  transmission  coefficient  T9  are  expected  to  coalesce  with  the  plane  wave  transmission 
coefficient  in  the  limit  of  an  infinitely-extended  (er  — ►  oo)  wavepacket.  The  distribution 
Pmax  as  a  function  of  energy  peaks  at  the  resonance  energy,  and  the  width  (FWHM)  of  this 
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Fig.  5.  Integrated  probability  P(t)  inside 
the  well  as  a  function  of  time.  The  mean 
energy  and  width  of  the  incident  packet 
are  the  energy  of  the  first  resonance  state 
and  400  A,  respectively.  The  dot-dashed 
line  is  an  extrapolation  of  P{t). 
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Fig.  7.  Pmax  in  the  upper  panel  and 
T(t  =  oo)  in  the  lower  panel  as  a  function 
of  the  mean  energy.  The  solid  line  shows 
the  plane  wave  transmission  coefficients. 
The  dashed,  dotted  and  dot-dashed  lines 
in  the  lower  panel  are  T 9  for  a  =  400,  500 
and  600  A,  respectively.  The  solid  circles, 
open  circles  and  solid  triangles  represent 
the  results  of  cr  =  400,  500  and  600  A, 
respectively. 
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Fig.  6.  Integrated  probability  P(t)  inside 
the  well  as  a  function  of  time.  The  mean 
energy  and  width  of  the  incident  packet 
axe  the  energy  of  the  second  resonance 
state  and  400  A,  respectively. 
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Fig.  8.  Build-up  time  as  a  function  of 
the  mean  energy.  The  explanation  of  the 
points  is  the  same  as  in  Fig.  7. 
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distribution  decreases  as  a  increases.  This  is  a  consequence  of  resonant  tunneling,  because 
for  the  case  of  on-resonance  scattering,  a  substantial  portion  of  the  wavepacket  tunnels 
through  the  first  barrier,  whereas  for  off-resonance,  most  of  the  packet  gets  reflected  by 
the  first  barrier. 

The  resonance  tunneling  process  can  further  be  demonstrated  by  looking  at  n  as 
a  function  of  packet  mean  energy,  as  displayed  in  Fig.  8  for  energies  in  the  vicinity 
of  the  second  resonance  state  and  a  =  400,  500  and  600  A.  This  tj  also  peaks  at  the 
resonance  energy,  which  can  not  be  explained  through  simple  dynamics  but  must  involve 
the  resonance  tunneling.  In  a  simple  dynamical  picture,  energetic  the  particle  takes  less 
time  to  to  travel  a  path  filled  with  obstacles  (repulsive  potential  profile  such  as  an  RTS), 
whereas  in  the  case  of  resonance  tunneling,  the  electron  with  energy  equal  to  the  resonance 
energy  bounces  back  and  forth  from  the  walls  of  the  well,  since  its  wavefunction  is  negligible 
at  the  walls  in  comparison  with  the  rest  of  the  places  in  the  well,  and  it  takes  much  more 
time  to  accumulate  probability  than  the  electrons  with  energies  off-resonance. 

In  conclusion,  we  have  shown  here  a  dynamical  picture  of  a  resonance  tunneling  pro¬ 
cess.  The  probability  inside  the  well  decays  exponentially,  and  excellent  agreement  is 
found  between  the  decay  constant  and  the  static  width  of  the  resonance.  We  have  also 
demonstrated  that  the  transmission  probability  coalesces  with  the  plane  wave  transmission 
coefficient  in  the  limit  of  an  infinitely-extended  wavepacket. 
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